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ABSTRACT 

In a previous paper (part I) , the mathematical properties of the cosmic microwave background radiation 
power spectrum which presents oscillations were discussed. Here, we discuss the physical interpretation: 
a power spectrum with oscillations is a rather normal characteristic expected from any fluid with clouds of 
overdensities that emits/absorb radiation or interact gravitationally with the photons, and with a finite 
range of sizes and distances for those clouds. The standard cosmological interpretation of "acoustic" 
peaks is just a particular case; peaks in the power spectrum might be generated in scenarios within some 
alternative cosmological model that have nothing to do with oscillations due to gravitational compression 
in a fluid. 

We also calculate the angular correlation function of the anisotropics from the WMAP-7yr and ACT 
data, in an attempt to derive the minimum number of parameters a polynomial function should have 
to fit it: a set of polynomial functions with a total of «6 free parameters, apart from the amplitude, is 
enough to reproduce the first two peaks. However, the standard model with six tunable free parameters 
also reproduces higher order peaks, giving the standard model a higher confidence. At present, while 
no simple function with six free parameters is found to give a fit as good as the one given by the 
standard cosmological model, we may consider the predictive power of the standard model beyond an 
instrumentalist approach (such as the Ptolemaic astronomy model of the orbits of the planets) . 
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1 Introduction 

Analysis of the Cosmic Microwave Background Radiation (CMBR) is one of the most important pieces 
of observational information supporting the standard model of cosmology, together with the abundance 
of light elements, the distribution of galaxies and their redshift-distance calibration, and other data. 
Moreover, in the last two decades, the CMBR has also become a source of accurate values for the 
parameters of this standard model of cosmology, which is usually called "precision Cosmology" . [1] The 
foundations of cosmology are thought to be definitively established, and now a quantitative science is 
pursued for which the fitting of the power spectrum of the distribution of CMBR anisotropics and few 
other data gives us the numerical values of the parameters in the equations governing the entire past, 
present and future of our Universe. 

The history of CMBR analysis is one of both successes and failures in the confirmation of predictions 
(that is, theoretical calculations carried out before the observations are available) and plenty of ad hoc 
fittings and corrections of the standard model. For instance, the fact that the CMBR presents an isotropic 
black body shape was a successful prediction. However, the first predictions of the (antenna) temperature 
of the CMBR were wrong: none of the predictions of the background temperature, which ranged between 
5 K and 50 K,[2] matched observations, the worst being the value of 50 K given by Gamow in 1961; in 
1965, the year of the discovery of the correct temperature, a value of 30 K was calculated by Dicke et 
al.[3] from the amount of helium production observed. The first predictions of CMBR anisotropics were 
also wrong, The first values predicted being one part in hundred or thousand; [4] however, this value could 
not fit the observations, which gave values hundreds of times smaller, so non-baryonic dark matter was 
introduced ad hoc to solve the question. 

CMBR analyses in the last decade have concentrated on anisotropies on small angular scales, smaller 
than the angular resolution of several degrees achieved by the COBE-DMR in the '90s. The power 
spectrum of the anisotropies of BOOMERANG and MAXIMA-1 showed a first peak at the Legendrc 
multipole t ~ 200, [5, 6] corresponding to angular scales of w 1.2°, which was interpreted as a discovery 
of the previously predicted acoustic peaks in a flat Universe with Q = Q m + Q\ = 1, with the greatest 
contribution from a dark energy component of fl\ = 0.7 (see review in Ref. [7] and references therein 
for further comments). Indeed, there are two parts in the comparison of theoretical predictions and 
observations: a successful one, and a half-successful/half- failure one. The totally successful prediction 
was the qualitative shape of the power spectrum containing several peaks: For instance, Peebles & Yu[8] 
had predicted a CMBR power spectrum with fluctuations, the acoustic peaks, from the hypothesis of 
primeval adiabatic perturbations in an expanding universe. However, the position of the first peak was 
not at the expected position when first measured: already in the mid-90s the position of the first peak 
was determined to be I « 200 with other experiments to measure small angular resolution anisotropies 
previous to BOOMERANG or MAXIMA-1 but with smaller sky coverage (Ref. [9] and references therein). 
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De Bernardis et al.[5] and Hanany et al.[6] just added a refinement in the measurement of the position of 
the first peak but they were not its discoverers. White et al.[9] realized that the preferred standard model 
of the time (an open Universe with il — ft m s=y 0.2 and without dark energy) did not fit the observations, 
so that they needed a larger il. This was one of the elements, together with SNIa observations, which 
would encourage cosmologists to include the new ad hoc clement: dark energy. Between 1997 and 2000 
this change of mentality in standard cosmology was produced, and then, in 2000, with the results of 
the new BOOMERANG and MAXIMA-1 experiments, cosmologists were proud to announce that new 
observations were giving exactly the results they expected. In any case, just paying attention to the 
first analyzes of the first peak in the mid-90s, we attribute a half-success to the prediction because, even 
though they failed to fit the observations with the preferred standard model of a curved open Universe at 
that time, the idea of a fiat Universe was also presented as a possibility in the '90s, and indeed a possibility 
preferred by inflationary paradigms. The position of the other peaks would also serve to constrain the 
cosmological models. The acoustic peaks on angular scales of 1 deg and 0.3 deg were predicted with the 
second peak nearly as high as the first one in lit + l)CV/(27r).[10, 11] Later, as data became available, 
its amplitude was reduced and the positions of the other peaks and their relative heights were also 
constrained from the model with a higher agreement with the data (e.g., Rcf. [12]). 

More recently, the analysis of CMBR polarization anisotropies (e.g., Ref. [13]) has also provided strong 
support for the standard cosmology. This indicates the way in which the CMBR is polarized due to the 
scattering of free electrons. Photon diffusion into regions of different temperatures are possible only when 
the plasma becomes sufficiently optically thin; these diffused photons could then scatter only while there 
are still free electrons left. Since photons could not diffuse too far, polarization cannot vary much over 
very large angular scales. 

Given this history of CMBR analyses, it is difficult to accept that our observations reproduce "predic- 
tions" . Nonetheless, whichever comes first, theory or observation, the fact remains that we have now have 
a cosmological model which is able to fit the CMBR power spectrum Ci quite accurately. It is usually 
thought that this may not happen by chance, given the apparent complexity of the power spectrum shape, 
whereas the six-parameter cosmological model is relatively simple. [13, 14] In order to fit the temperature- 
polarization cross power spectrum and the polarization-polarization power spectrum, [13] one would need 
an extra parameter (optical depth), so a total of at least seven free parameters are necessary. In this 
paper, we will not talk about the polarization anisotropies and their fits with a seven-parameter model. 
Roughly speaking, the relationship between temperature-temperature and the polarization-temperature, 
or polarization-polarization power spectra is expected since they arc different ways of seeing the same 
light with different filters; a full discussion of the generation of polarization is not the subject of this 
article. 

Some critics of the standard model [15, 16] claim that there is little value in the fact that the standard 
cosmological model might fit some observational cosmological data if the number of free parameters in 
the model is comparable to the number of independent parameters characterizing the observations. For 
instance, Ptolemaic geocentric astronomy may fit the observations of the orbits of the planets but with 
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too high a number of free parameters in the theory. The principle of Occam's Razor tells us that a theory 
is better when the number of free parameters is low. Occam's Razor can indeed be understood with a 
Bayesian analysis, in which it is tested how probable a model is with a given number of free parameters 
with respect to other models with a higher number of free parameters. [17, 18] Otherwise, if a theory 
has a high number of free parameters, it loses credibility because it is always possible to create a "false" 
model to fit some data when the number of free parameters is comparable to the number of degrees of 
freedom in the data. Philosophers of science, when talking about cosmology, associate this approach 
with "instrumcntalism" (e.g., Ref. [19]). And saying that the power spectrum contains hundreds of 
independent parameters for a given resolution is not correct, because the different values of Cg for each 
t are not independent in the same sense that hundreds of observations of the position and velocity of a 
planet do not indicate hundreds of independent parameters, the information on the orbit of planet being 
instead reduced to only six Keplerian parameters. 

For the reasons just given, understanding how much information is in the power spectrum is important 
for key questions in the discussion of the fundamentals of cosmology, and this is the matter we shall discuss 
in this paper. In spite of the huge number of papers dedicated to analyzing how to fit the power spectrum 
with standard cosmological models, there is a dearth of literature concerning the properties of this power 
spectrum independently of the standard physical model. Indeed, there are here two main points which 
should be clarified: 

1. Are the oscillations something atypical in a power spectrum? That is, should we consider the fact 
the power spectrum contains oscillations a successful prediction of the standard cosmological model 
which cannot be produced by any other means? 

2. How much information is contained in the power spectrum? That is, how many free parameters in 
a function are necessary to fit the CMBR power spectrum? Should we consider the fitting of the 
power spectrum with a model of six free parameters as a validation of the standard cosmological 
model? 

Our paper is structured as follows. In §2, we discuss the physical origin of the oscillating peaks in 
the power spectrum, and of the corresponding scales in the angular correlation function, which corre- 
sponds to the physical interpretation of the mathematical toy models developed by Lopez-Corredoira & 
Gabrielli [20] (hereafter Paper I). In §3, we analyze the angular correlation function of CMBR anisotropics 
on small scales, and we find that a set of polynomial functions with a total of «6 free parameters may 
roughly fit the angular correlation function. Finally, in §4 we will discuss these results and summarize 
them in §5. 
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2 The physical origin of oscillations in the power spectrum 

The fluctuation field on the sky, %f(6, 4>) can be decomposed into spherical harmonics in the sphere: 

s-rp oo +f- 

e=o m=-l 

Assuming statistical isotropy, the angular power spectrum Ce is defined as 

C t = (|a, m | 2 ) = ^ £ l a ^| 2 ' ( 2 ) 

m=—£ 

whereas the angular two-point correlation function is defined as (Paper I): 

(3) 

where 6 = 2 — 0i. The two-point correlation function can be expanded in terms of Legendre polynomials 

1 oo 

C W = T + l)C t P t [coB(0)] , (4) 

where Pi are the Legendre polynomials. By inverting Eq. (4) we find[21] 

G\ = 2ir f C(9)P l [cos(9)]sm(9)d9. (5) 
Jo 

A remarkable advantage of the use of self-correlation rather than the power spectrum is that it gives 
more direct information on the physical space (Paper I), so we shall use it here. 

2.1 CMBR anisotropies in terms of the Toy Models of Paper I 

In Paper I, it was shown that oscillations observed in Ce could be generated by a self-correlation C(9) 
with an abrupt transition for some angle. That is, for an angular correlation function such that 

f gi(0), 9<9 } 
C{9) ' " ) , (6) 

( g 2 (9), 9 > 9 J 

with <7i"^(#o~) 7^ g^i^o) an d n somc integer number greater than or equal to one; that is, a function 
with some derivative which is not continuous at 9o- Also in Paper I, we saw that a toy model in which we 
distribute disks of fixed or variable radii with some distribution of temperature T, but with the value of 
the field outside the disks being T = 0, we get angular correlation of the type of Eq. (6) and consequently 
we get oscillations in the power spectrum. We wonder now whether these toy models of Paper I could 
produce a CMBR distribution which mimicked the observed power spectrum/angular correlation, and 
what their physical interpretation would be. 

We could interpret the disks as projections of spherical clouds in 3D space that emit/absorb radiation 
or interact gravitationally with it. Toy models 1 of Paper I restrict the angular radius of the disks/clouds 
to be a constant. The case of toy model lb, with the restriction that disks do not intersect in their 
projection, produces a dip for 2R < 9 < 2R + S for some small S. Other examples of matter distribution 
over spherical bubbles that do not produce oscillations are also possible (e.g., Ref. [22], §13.4). 
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Just looking at Fig. 1/top, we can see there are two clear regimes for C(9) in the real data with 
a transition at 9 w 1.2°, which is directly comparable with the toy models of Paper I. A value of 
R w 0.6° would be necessary to assimilate it to the results of the toy model 1 (Fig. 3 of Paper I). 
The case of toy model lb, with the restriction that disks do not intersect in their projection, produces 
a dip for 2R < 9 < 2R + 5 for some small 6, and this is indeed observed in the real CMBR self- 
correlation too (see Fig. 1/Middle). From this analysis, we get the idea that with a simple model of a 
distribution of clouds of constant physical radius at a fixed distance, with some anticorrelation among 
them, which emit / absorb radiation or interact gravitationally with the radiation, we can get the CMBR 
self-correlation, its corresponding power spectrum with the oscillations, or even with a distribution of 
clouds of variable radii and distances (toy model 2). 

The standard model obeys the scenario of the toy model 1 of Paper I: The redshift associated with 
the matter-radiation decoupling gives us a fixed distance and the constant physical size is given by the 
size of the acoustic horizon. The origin of oscillations in the power spectrum is well understood in terms 
of a photon-baryon fluid compressed by gravitational attraction produced by local density fluctuations 
giving rise to acoustic oscillations. [8], [23] (§6. 6. 2) The angular size 9 = 1.2° corresponds to the acoustic 
horizon size at recombination. [24, 25] In any case, apart from the standard model, any physical scenario 
resembling the conditions of the toy models 1 may produce a power spectrum with the observed features. 

If we were to think about non-standard cosmologies, considering a distribution of disks with a contant 
radius at constant distance would be very restrictive. But toy models 2 of Paper I teach us that a 
constant radius and distance are indeed unnecessary for those clouds to generate an abrupt transition. 
The case of toy model 2a shows us that a distribution of variable sizes of the clouds still produces abrupt 
transitions at the minimum (2R m i n ) and maximum (2R max ) angular sizes of the cloud. The case of toy 
model 2b shows us that, for a fixed linear size (L), and a distribution of clouds in 3D space in a range of 
distances between r m i n and r maxi we get abrupt transitions in the correlation function at 2L/r max and 
2L/r m i n . A combination of variable sizes and variable distance is also possible with similar conclusions. 
A distribution of clouds over an infinite distance, apart from the problem of Olber's paradox, would 
produce a value of C{9) that is continuous in all of its derivatives and that would not give rise to any 
oscillation. 

The conclusion to take from the analysis of these toy models, therefore, is that we could in principle 
generate a distribution of CMBR anisotropies with oscillations in the power spectrum by setting a fluid 
with some clouds of overdensities (with respect to the average density) , which emit / absorb radiation or 
interact gravitationally with the radiation over a finite range of sizes and distances. 

The question which now arises is how much fine tuning we need to fit a model following the mathemat- 
ical recipes of the toy models in Paper I onto the observed CMBR anisotropy distribution. Indeed, the 
toy models of Paper I have some tunable parameters and functions f(9), u>(6) or A(R), whose variation 
can generate the different C(0) and Cg with oscillations: 




(7) 
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h(0) = h[f(6),w(6)}, 

where h(9) is the correlation function derived from toy model 1 for disks of angular radius i?o, and its 
integration for different angular radii with amplitudes A(R) gives rise to toy model 2. The function f(9) 
gives the isotropic distribution of antenna temperature within each disk. The function uj(9) gives the 
two-point correlation function in the distribution of the centers of the disks. 

One could produce a distribution of disks as in toy models 1 and 2 to produce such fits. Theoretically, 
a fit of C{9) through the toy-model representation would get the best functions f(9), uj(9) or A(R), or 
the different families of functions if there is degeneracy. From Eq. (7), it looks generally possible to 
recover functions A, /, uj, which can produce an observed C{9). Moreover, if we fix a toy model 1 h(9), 
then, in principle, we could find a function A(R) to fit the data by inverting Eq. (7). However, in such 
a case this function A{R) might be very complex and difficult to parameterize with few parameters, so 
a free combination of the three functions should be allowed to get the fit with a minimum number of 
parameters. 

Nonetheless, in practice this fit cannot be carried out directly. Testing all the different combinations 
of functions would be a very time-consuming computer calculation: for each set of /, uj, the calculation 
of h takes ~10 minutes for 9 < 5° in normal computers with the resolution given in Paper I, and if we go 
to much higher 9 (required for <C 1; the computing time for a given 9 is larger for larger 9) the time is 
increased to hours for each h, so exploring the multiple possible combinations of / and uj shape functions 
would require an excessive time, even with a super-computer. Instead, for the question of the number 
of necessary parameters, we will analyze in the following section the fit of C{9) by a generic polynomial 
function. 

3 Parameterization of CMBR angular correlation function 

Our problem is to determine the amount of independent information contained in C(9), or its Fourier 
transform P(k) or C;. The functional dependence is simpler for the real-space, C(9), than in frequency 
space (paper I), so we will concentrate on the analysis of the former. 

In signal processing of multidimensional signals, for example in computer vision, the intrinsic dimension 
of the signal describes how many variables are needed to represent the signal. For a signal of N variables, 
its intrinsic dimension is an integer number M < N . The intrinsic topological dimensionality of data 
refers to the minimum number of free parameters needed to generate a data set. [26] This concept of 
"intrinsic dimension" can also be applied to our problem: knowing the function with a minimum number 
of free parameters which fits our C{9). We want to know the intrinsic dimension of C{9). It contains 
hundreds of values for different 9, for the given resolution, but they are not independent. However, 
although the problem is defined in a precise way, the ways of finding the solution are not. There is some 
literature on algorithms to get the intrinsic dimension on data structures in high-dimensional space (e.g., 
Refs. [27, 28]), but these are not applicable to our specific problem. 

In this paper, we use fits to simple polynomial functions as a way to estimate a number of free 
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parameters, but the real intrinsic dimension of C(0) might be lower with other more complex functions. 
Certainly, there may exist another non-trivial function shape which produces better fits. But we prefer 
not to explore all the possible known functional shapes because this would reduce the probability to being 
a chance fit. That is, the larger the number of test functions that we use, the higher the probability of 
obtaining a better fit, but this probability should be divided by the number of total test functions used 
in the trials. 

3.1 WMAP angular correlation function 

The direct valculation of the angular correlation function from CMBR data is subject to several difficulties. [25] 
Rather, we will carry out the calculation of the self-correlation of WMAP (Wilkinson Microwave Anisotropy 
Probe)-7yr on small scales corresponding to the transform [Eq. (4)] of the published power spectrum. [13] 
The result is shown in Fig. 1, which illustrates the combined power spectrum using all the available 
frequencies. Indeed, when we talk about data from observations, C° bs - = WeC]! eal , with Wg as a window 
function, usually with a shape We = exp[—£(£+l)cr 2 } 7 with a related to the size of the beam. Throughout 
this paper we shall concentrate on the "real" power spectrum, once corrected for the window function, so 
from now on Cg = C]i eal , and the correlation function C(0) will also refer to the real correlation function 
rather than that directly obtained from observations. We also assume that foreground contamination has 
been correctly subtracted, regardless of doubts that may be entertained concerning this subtraction. [29] 
The error bars stem from the error in the power spectrum due to the measurement (error derived from 
diagonal terms of Fisher matrix; the multipole moments are slightly coupled, so a correct treatment of 
errors would require use of the entire Fisher matrix; small, correlated contribution due to beam and 
point source subtraction uncertainty are not included; cosmic variance error is not included), calculated 
through Monte Carlo simulations. The total errors are expected to be slightly larger but not much, as 
we will see below. 

We would like to determine how many free parameters are necessary to describe C(0) in Fig. 1. In 
Fig. 1, there are best weighted fits of this correlation function within 0.2° < < 5° with two functions of 
class C 1 (the function and its first derivative are continuous) respectively with 4 and 6 free parameters 
(apart from the amplitude): f± and fe, which are respectively consecutive polynomials of second degree 
for two or three different angular ranges satisfying the continuity of fi and /-. 

f e [9(deg.)}^K 2 ) = (8) 

{ 12067 - 161180 + 67460 2 , 9 < 1.239 
2452 + 593.401 - 35420?, 1.239 < 6 < 1.391 >, 
2460 — 485.16*2 + 53.240f, 9 > 1.391 
6>i = 0- 1.239; 2 = 9 - 1.391 

[Fit with N = 192 points giving \ 2 = 584]. 

Note that all of the numbers given in Eq. (8) are not independent. In Eq. (8), three numbers for each 
of the three second-degree polynomials are given, and two angles where the transitions are produced, 
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Figure 1: Top/Middle: Self-correlation function of WMAP-7yr CMBR angular anisotropies and fits to 
it. Errors of WMAP correlation only include instrumental noise. Bottom: Power spectrum of these fits, 
together with the WMAP-7yr data. Also included is the best fit with the ACDM mode and its angular 
correlation (plus a constant term — 401/iK 2 in order to compensate for the effects of cosmic variance at 
low i in the data). 
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so there are eleven numbers. But imposing continuity on the function and its first derivative, reduces 
the number of independent parameters to four (two for each transition point). The amplitude was not 
included among the free parameters, so in total we have six free parameters in fa. 

U[9(deg.)](pK 2 ) = (9) 

12144- 164340 + 7O120 2 , 9 < 1.140 1 
2523 - 446.76*3 + 39.290§, 9 > 1.140 J ' 
9 3 = 9- 1.140 

[Fit with N = 192 points giving X 2 = 1844]. 

Another function of class C° (continuity in the function but not in its derivative) with four free 
parameters which roughly fits C(9) is: 

94 [9(deg.)](nK 2 ) = (10) 

12073 - 161440 + 67590 2 , 9 < 1.320 1 

2539-877.4[Zn(0)-Zn(1.32O)], 9 > 1.320 J 
[Fit with N = 192 points giving X 2 = 768]. 

We plot these fits in Fig. 1/top- middle. We also plot the fit given by the standard A-CDM model, 
plus a constant term Ck in order to compensate for the effects of cosmic variance at low I in the data. 
It is this term that minimizes the chi-square in the fit: Ck = —401 ^K 2 , \ 2 = 330 in the same N = 192 
points. The fact that we obtain X 2 /N ~ 1.7 instead of one with the standard model indicates that the 
error bars may be underestimated by «25%. 

The fits with fa and g 4 are better than that with fa; they fit the small dip in the transition at 9 rts 1.2° 
better; in any case, fa also gives a fairly acceptable fit. Bashinsky & Bertshinger[24] showed that, if this 
dip is cut out, the oscillations in the power spectrum after the first peak would be removed; but, as we 
have shown, this depends on how the dip is cut out. Other minor oscillations of the correlation function 
are most probably fluctuations due to noise. 

The power spectrum associated with the fitted correlations fa, fa and g\ (for fa and fa, over 5 degrees, 
we have extrapolated linearly till the value of the function reaches null value, and zero beyond), derived 
through Eq. (5) applied to C — fa, C — fa, C — 54 respectively, is plotted in Fig. 1/bottom. The 
agreement with the observed power spectrum in the first two peaks is acceptable. All of them reproduce 
the first peak, and the second peak is more or less well fitted for fa or g±. The amplitude of the third 
peak, however, is not well fitted, except for its center and width. 

The minor oscillations in the power spectra of Fig. 1/bottom are produced by the method of extrap- 
olation. We can perform different extrapolations of the same fa for 5° < 9 < 0*, all of them main- 
taining a continuous function and derivative in 9 = 5°, such that the extrapolation reaches fa{6*) = 
and we set fa(9 > 0*) = 0: 1) Extrapolation 1 (the linear one used for Fig. 1/down): fa[0(°) > 
5] = /e(5) + #(5) x (9 - 5); 2) Extrapolation 2: / 6 [0(°) > 5] = / 6 (5) + |/£(5) x [(9 - 4) 2 - 1]; 3) 
Extrapolation 3: fa[9{°) > 5] = / 6 (5) +2/^(5) x [y/9~ 4 - 1]; 4) Extrapolation 4: fa[9{°) > 5] = 
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Figure 2: Power spectrum (Fourier transform) of f§ with different extrapolations for 5° < 9 < 9* such 
that the extrapolation reaches fe(0*) = and we set fe(9 > 6*) = 0: 1) Extrapolation 1 (the linear 
one used for Fig. 1/down): f 6 [6(°) > 5] = / 6 (5) + /^(5) x (9 - 5); 2) Extrapolation 2: f 6 [6(°) > 5] = 
/e(5) + yi(5) x [(9 - 4) 2 - 1]; 3) Extrapolation 3: f 6 [6(°) > 5] = / 6 (5) + 2/^(5) x [^9~ 4 - 1]; 4) 
Extrapolation 4: / 6 [0(°) > 5] = / 6 (5) + 2000 x [exp [5^(5) x {9 - 5)] - l] . 

fe(5) + 2000 x [exp [^c^Jo/^^) x (9 — 5)] — l] . Fig. 2 shows these four cases. As can be seen, the re- 
sult of the inversion is the same, and the change of extrapolation only affects to the oscillations of high 
frequency. 

Although the power spectrum seems to contain a lot of information with all their oscillations, indeed 
all the information is concentrated in the abrupt transition of the self-correlation function at scales of 
1.1-1.4 degrees. In our examples, the second derivatives of fe are not continuous at 9 — 1.239° and 
9 = 1.391°, the second derivative of f± is not continuous at 9 = 1.140°, and the first derivative of 54 is 
not continuous at 9 — 1.320°. A simple function with six free parameters is able to reproduce more or 
less the correlation. Even with four free parameters, a rough fit is obtained, especially for 34. The best 
fit with the standard cosmological model[13] (see Fig. 1/down) is indubitably better. The last WMAP 
data of the 7-yr release have precisely allowed a better determination of this third peak, [13] being almost 
with the same amplitude as the second peak, and this gives an important advantage to the standard 
cosmology with respect to a generic polynomial fit, which tends to reduce the power of the peaks in 
their consecutive sequence. Hu et al.[12] also showed with BOOMERANG and MAXIMA-1 data that 
the height of the second and third peaks are very similar. The best value of x 2 /N for our fits was 3.0 for 
fs, which indicates that the standard model (x 2 /N = 1.7) produces a better fit. 

The fact that we have found a polynomial fit such as fe with six parameters giving an approximate 
fit does not mean that this is the best possible fit with generic functions with < 6 free parameters. The 
intrinsic dimension, might be 6 or even lower and get a much better fit with other kinds of functions. 
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Figure 3: Top: Self-correlation function of WMAP-7yr+ACT CMBR angular anisotropics and fits. Bot- 
tom: Power spectrum of this fit, together with the WMAP-7yr and ACT data. 

3.2 WMAP+ACT angular correlation function 

Higher-order peaks were obtained in other experiments, for instance in the ACT (Atacama Cosmol- 
ogy Telescope [30]) project, and these data might be used for an analysis like the one carried out in 

the previous subsection. Here, I will carry out the analysis with the 148 GHz data once corrected 

c obs - 

for the window function, that is: C\ = ^ for 500 < I < 9750 (using the available points at 
http:/ /lambda. gsfc.nasa.gov/product/act/act_spectra_get.cfm and with linear interpolation for the rest); 
for I < 500 we will use the WMAP data given in the previous subsection. Figure 3 plots the self- 
correlation of this WMAP+ACT power spectrum at small scales, corresponding to its transform given 
by Eq. (4). We also use Monte Carlo simulations to estimate the error bars from the errors given by the 
analytical expression from Rcf. [30] . The results are very similar to the WMAP self-correlation. 

Again, a best-weighted fit of this correlation function is carried out, within 0.2° < 9 < 5° with a 
function of class C 1 and six free parameters, which is a consecutive polynomial of second degree for three 
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different angular ranges satisfying the continuity of /j and f[: 

f e , A [0(deg.)}(^K 2 ) = (11) 

{ 12159- 163490 + 68780 2 , 6 < 1.213 
2448 + 335.66*1 - 16720?, 1.213 < 6 < 1.457 >, 
2430 - 481.902 + 54.1202, > 1.457 

0i ee 0- 1.213; 2 = - 1.457 

[Fit with N = 192 points giving X 2 = 260]. 

This fit is much better than the equivalent one with WMAP data only, with x 2 /N = 1.35. Nonetheless, 
a new transform of this correlation function to obtain the power spectrum (see Fig. 3/below) has the 
same problem as with the WMAP data only. The third and subsequent peaks are more or less fitted in 
position but not in amplitude. Again, the standard model gives a much better fit. [30] Note, however, 
the ACT works at 148 and 218 GHz, from which a fitting model needs to include lensed CMB, Sunyaev- 
Zel'dovich, and foreground contributions with a total of ten or more free parameters. The intergalactic 
medium might become dominant on I >^ 2000. [31] 



4 Discussion 

Apart from the standard model, any physical scenario resembling the conditions of the toy models 1 or 2 of 
Paper I, with some sets of parameters and functions /(0), w(0) or A(R) which fit the correlation of Fig. 1, 
may produce a power spectrum with the observed features. Nevertheless, we must bear in mind that any 
possible alternative scenario should also explain the almost perfect black body shape of that radiation, 
and the value of around 3 K in its antenna temperature, for which there are exotic ideas in the literature, 
but which are not free from caveats and inconsistencies. There are, for instance, proposals in terms of a 
local origin of CMBR and its thermalization of stellar radiation by the intergalactic medium in the form 
of thermalization produced by an aether, [32] carbono needles ( "whiskers" [33, 34]), cosmic meteoroids 
("cosmoids" [35]), radiation of electrons orbiting protons in the atoms [against quantum mechanics] , [36] 
or electrons or ions upon the action of magnetic fields. [37, 38, 39, 40, 41] The origin must be very local 
in order to preserve the black-body shape, or to be non-redshifted, which seems problematic. Some 
proposals have even placed the origin of the CMBR within our own Galaxy, in the local bubble within 
100 pc,[42] but it is puzzling to understand by this means the Sunyaev-Zel'dovich imprint of clusters of 
galaxies in the CMBR [43] or any possible correlation of the CMBR with galaxies. 1 

Another relevant feature of CMBR anisotropics is the gaussianity of its fluctuations. However, this 
should not be a major problem for alternative cosmologies since many different origins of the fluctuations 
are expected to be Gaussian, although not all of them. Theories involving inflation generally predict a 
pattern of Gaussian noise, whereas theories based on symmetry breaking and the generation of defects 

1 This has indeed not been found yet, since the measurements of the cross-correlation of WMAP with galaxy surveys are 
not significant; see Ref. [44] and references therein. 
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have more distinctive signatures. [45] Moreover, it is not totally clear yet whether the CMBR fluctuations 
are exactly Gaussian: A non-Gaussian distribution has been claimed by several authors (e.g., Ref. [29], 
§2.3 and references therein; Refs. [46, 47, 48, 49]); however, other analyses[50, 51] claim that only 
a few regions have such non-Gaussian anisotropies due to contamination, e.g. the Corona Borealis 
supercluster region, and most of the regions in the sky are Gaussian. It has also been argucd[52, 53] that 
the non-gaussianity is associated with cold spots of unsubtracted foregrounds; even the lowest spherical 
harmonic modes, which should be the cleanest, in the map are significantly contaminated with foreground 
radiation. [54] 

As said in §3, a simple function with six free parameters is able more or less to reproduce the first two 
peaks of the power spectrum. Even with four free-parameters, 54, a reasonable fit is obtained. Hence, 
one would expect that any theory of the type of toy models 1 or 2 with this number of free parame- 
ters could represent the data. The standard cosmological model fits the power spectrum with six free 
parameters. 2 [14, 13] This is astonishing but it would be much more surprising if it could fit it without 
free parameters or with only one or two free parameters fitted to the power spectrum. Certainly, the 
same parameters which are used to fit the CMBR power spectrum are able to fit other cosmological 
data, but more than six parameters are necessary, as well as additional information such as initial condi- 
tions, conditions of stellar formation, galaxy formation, how dark matter distributes in galaxies, etc. A 
global analysis of the cosmological models would require the examination of all of the available indepen- 
dent sources of cosmological data (nucleosynthesis, supernovae, gravitational lensing, etc.), and to check 
whether they are comparable to the number of free parameters in the model. We do not carry out this 
analysis in the present paper. At present, we refer only to the CMBR anisotropies. 

Observe that we are not claiming here that fe or any best fit could substitute the standard theory in 
order to explain the observed Ci or C{6). The function j% is just a simple fit which does not represent 
any physical model. The toy models of Paper I may represent a physical model, but not the polynomial 
functions of §3. A simple fitting recipe cannot substitute a physical theory even if it requires a smaller 
number of parameters to reproduce experimental/observational data. A physical theory can be substi- 
tuted only by another theory which is able to explain all the phenomena taken into account by the first 
one and requires fewer subtleties and ad hoc assumptions. This is nothing else that the application of 
the principle of Occam's razor, but this is valid for homogeneous things, i.e. two theories (or physical 
explanations) and two fitting procedures (for instance the comparison between a fit with a sum of expo- 
nentials or a sum of power laws) and not a theory and a fitting procedure that have different intrinsic 
value. Nonetheless, our aim here was just to show the capacity that a wrong theory would have to fit 
C{&) with six free parameters by chance, and we have shown that up to now [further research is necessary 
in the fitting of C{9)\ the standard cosmological model makes it better than would be expected if it were 
a wrong theory. That is, the standard model looks correct, otherwise it would not be able to produce the 

2 Plus many other parameters which introduce second-order changes. And, even so, there is a degeneracy in the solutions 
with different values of -Ho an d CI a.: WMAP data and the large scale structure of galaxies could be reproduced without 
explicitly requesting the existence of dark energy, [55] i.e. with A = 0. This degeneracy is broken by adding cosmological 
information from other sources, for instance, from SNIa data. 
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accurate fit of C(6). If wc had found that produced a fit as good as or better than the standard model 
then it would not mean that a fit f§ represents a theory which could be an alternative to the standard 
model, but we might wonder whether standard cosmology is real or just an instrumental model to save 
the data (as in Ptolemaic astronomy), at least with regard to the CMBR anisotropics. However, this was 
not the case. 

In cosmologies different from the standard model, there were also attempts to fit the CMB power 
spectrum with models with a similar number of free parameters. Narlikar et al.[56, 57] fitted the Ci with 
five-six parameters apart from the amplitude within a model which has nothing to do with the origin of 
the CMBR in the standard cosmology. Power spectrum peaks at I w 6 to 10, 180 to 220, and 600 to 900 
are shown to be respectively related in their Quasi-Steady State cosmology to curvature effects at the 
last minimum of the scale factor, clusters, and groups of galaxies. For a MOND (Modified Newtonian 
Dynamics) cosmology without cold dark matter, it is also possible with few parameters to fit the power 
spectrum. [58, 59] The standard model still fits the third peak better than our generic polynomial fit or 
most alternative cosmological scenarios as a function of the same number of free parameters, although 
there are also alternative scenarios (e.g., Ref. [59]) in which the fit is equally good. 

5 Conclusions 

Summing up, any hypothesis which tries to explain the origin of CMBR anisotropies must explain the 
oscillations of the power spectrum, apart — of course — from the black-body emission. Here and in Paper 
I, it has been shown that it is a rather normal characteristic expected from any fluid with clouds of 
overdensities that emit/absorb radiation or interact gravitationally with the photons, and with a finite 
range of sizes and distances for those clouds. Apart from the standard cosmological model, other scenarios 
may also follow these conditions. The interpretation of "acoustic" peaks is just a particular case; peaks 
in the power spectrum may be generated in scenarios which have nothing to do with oscillations due to 
gravitational compression. 

The angular correlation function is well fitted by a set of polynomial functions with a total of «6 
free parameters, apart from the amplitude. In the power spectrum, the first two peaks are more or less 
reproduced with these fits, whereas the third peak and higher-order peaks are not. In this respect, the 
standard cosmological model, or the Quasi-Steady State Model (both of them using six free parameters 
to fit the CMBR power spectrum) have not got too much predictive power to describe the first two peaks 
properties. Nevertheless, the last WMAP data or ACT data have given us the chance to test the standard 
theory with a good fit of the third peak too, providing greater confidence in the standard model. 

Further research is necessary along these lines in order to sec whether fits of the CMBR correlation 
function/power spectrum as good as the standard model one can be produced with totally different 
models and a similar number of free parameters. 

Acknowledgments: 



15 



Thanks are given to F. Sylos Labini, A. Gabrielli, F. Atrio-Barandela and the anonymous referee for 
helpful comments. Thanks are given to T. J. Mahoney for proof-reading this paper. WMAP is the result 
of a partnership between Princeton University and NASA's Goddard Space Flight Center. ACT operates 
in the Chajnantor Science Preserve in northern Chile under the auspices of the Comision Nacional de 
Investigation Cienti'fica y Tecnologica (CONFCYT). 

References 

[1] J. R. Primack, New Astr. Rev. 49 (2005) 25. 

[2] T. Van Flandern, Reflector (Astronomical League) XLV (1993) 4. 

[3] R. H. Dicke et at, Astrophys. J. 142 (1965) 414. 

[4] R. K. Sachs and A. M. Wolfe, Astrophys. J. 147 (1967) 73. 

[5] P. de Bernardis et al, Nature 404 (2000) 955. 

[6] S. Hanany et al, Astrophys. J. 545 (2000) L5. 

[7] W. Hu and S. Dodelson, Ann. Rev. Astron. Astrophys. 40 (2002) 171. 

[8] P. J. E. Peebles and J. T. Yu, Astrophys. J. 162 (1970) 815. 

[9] M. White et al, Mon. Not. R. Astron. Soc. 283 (2006) 107. 

[10] J. R. Bond and G. Efstathiou, Mon. Not. R. Astron. Soc. 226 (1987) 655. 

[11] H. E. Jorgenscn et al, Astron. Astrophys. 294 (1995) 639. 

[12] W. Hu et al, Astrophys. J. 549 (2001) 669. 

[13] D. Larson et al, Astrophys. J. Suppl. Ser. 192 (2011) 16. 

[14] H. Pciris, Contemporary Phys. 46(2) (2005) 77. 

[15] M. J. Disney, American Scientist 95(5) (2007) 383. 

[16] M. Lopez-Corredoira, Sociology of Modern Cosmology, in Cosmology Across Cultures (ASP Conf. 
Ser. 409), eds. J. A. Rubino-Martfn et al. (Astron. Soc. Pacific, S. Francisco), p. 66. 

[17] H. Jefferys and J. Bcrger, American Scientist 80(1) (1992) 64. 

[18] J. O. Bcrger and W. H. Jefferys, J. Italian Statistical Soc. 1 (1992) 17. 

[19] F. J. Solcr Gil, Themata 40 (2008) 117. 

[20] M. Lopcz-Corrcdoira and A. Gabrielli, Physica A 392 (2013) 474 (Paper I). 

16 



[21] C. J. Copi et al, Advances in Astronomy, id. 847541 (2010). 

[22] P. Coles and F. Lucchin, Cosmology. The Origin and Evolution of Cosmic Structure (John Wiley & 
Sons, London, 2002). 

[23] A. Gabriclli et al., Statistical Physics for Cosmic Structures (Springer, Berlin, 2005). 
[24] S. Bashinsky and E. Bertshinger, Phys. Rev. Lett. 87(8) (2002) 1301. 

[25] C. Hernandez-Monteagudo, A. Kashlinsky and F. Atrio-Barandela, Astron. Astrophys. 413 (2004) 
833. 

[26] A. K. Jain and R. C. Dubes, Algorithms for clustering data (Prentice Hall, Englewood Cliffs, 1988). 

[27] B. Kegl, Neural Information Processing Systems 16 (2002) 681. 

[28] T. Ban, C. Zhang and S. Abe, Journal of Intelligent Systems 33(1) (2008) 3. 

[29] M. Lopez-Corredoira, J. Astrophys. Astron. 28 (2007) 101. 

[30] S. Das et at, Astrophys. J. 729 (2011) 62. 

[31] F. Atrio-Barandela and J. P. Miicket, Astrophys. J. 643 (2006) 1. 

[32] S. V. M. Clubc, Mon. Not. R. Astron. Soc. 193 (1980) 385. 

[33] F. Hoyle, G. Burbidge and J. V. Narlikar, Astrophys. J. 410 (1993) 437. 

[34] F. Hoyle, G. Burbidge and J. V. Narlikar, Mon. Not. R. Astron. Soc. 267 (1994) 1007. 

[35] R. K. Soberman and M. Dubin, Dark Matter is Baryons [arXiv:astro-ph/0107550]. 

[36] G. P. Shpenkov and G. Kreidik, Revista Ciencias Exatas e Naturais 4 (2002) 9. 

[37] E. J. Lerner, Laser and Particle Beams 6 (1988) 457. 

[38] E. J. Lerner, Astrophys. Space Sci. 227 (1995) 61. 

[39] A. Brynjolfsson, Redshift of photons penetrating a hot plasma [arXiv:astro-ph/0401420]. 

[40] V. Krishan, Optical Depth of the Cosmic Microwave Background due to Scattering and Absorption 
[arXiv:0909.0125]. 

[41] D. Crawford, J. Cosmology 13 (2011) 3947. 

[42] J. V. Narlikar, J.-C. Pecker and N. C. Wickramasinghe, The Local Contribution to the Microwave 
Background Radiation(MBR), in Evolution of Cosmic Objects through their Physical Activity, eds. 
H. Harutyunian, A. Mickaelian and Y. Terzian (Gitutyun Publ. House of NAS RA, Yerevan), p. 259. 

[43] N. Aghanim et al, Astron. Astrophys. 536 (2011) A10. 



17 



[44] M. Lopez-Corredoira, F. Sylos-Labini and J. Betancort-Rijo, Astron. Astrophys. 513 (2010) A3. 
[45] D. Coulson et al, Nature 368 (1994) 27. 

[46] C. Racth, P. Schuecker and A. J. Banday, Mon. Not. R. Astron. Soc. 380 (2007) 466. 

[47] A. Bernui, C. Tsallis and T. Villela, Europhysics Letters 78(1) (2007) 19001. 

[48] J. D. McEwen et al, Mon. Not. R. Astron. Soc. 388 (2008) 659. 

[49] G. Rossi et al, Mon. Not. R. Astron. Soc. 399 (2009) 304. 

[50] J. A. Rubino-Martm et al, Mon. Not. R. Astron. Soc. 369 (2006) 909. 

[51] J. D. McEwen et al, Mon. Not. R. Astron. Soc. 371 (2006) L50. 

[52] X. Liu and S. N. Zhang, Astrophys. J. 633 (2005) 542. 

[53] R. Tojeiro et al, Mon. Not. R. Astron. Soc. 365 (2006) 265. 

[54] L.-Y. Chiang, P. D. Naselsky and P. Coles, Astrophys. J. 664 (2007) 8. 

[55] A. Blanchard et al, Astron. Astrophys. 412 (2003) 35. 

[56] J. V. Narlikar et al, Astrophys. J. 585 (2003) 1. 

[57] J. V. Narlikar, G. Burbidge and R. G. Vishwakarma, J. Astrophys. Astron. 28 (2007) 67. 
[58] S. S. McGaugh, Astrophys. J. 611 (2004) 26. 

[59] G. W. Angus and A. Diaferio, Mon. Not. R. Astron. Soc. 417 (2011), 941. 



18 



